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Abstract 

Dynamical symmetries of Hamiltonians quantized models of dis- 
crete non-linear Schrodinger chain (DNLS) and of Ablowitz-Ladik 
chain (AL) are studied. It is shown that for n-sites the dynamical 
algebra of DNLS Hamilton operator is given by the su(n) algebra, 
while the respective symmetry for the AL case is the quantum al- 
gebra su q (n). The g-deformation of the dynamical symmetry in the 
AL model is due to the non-canonical oscillator-like structure of the 
raising and lowering operators at each site. 

Invariants of motions are found in terms of Casimir central ele- 
ments of su(n) and su q (n) algebra generators, for the DNLS and QAL 
cases respectively. Utilizing the representation theory of the symmetry 
algebras we specialize to the n = 2 quantum dimer case and formulate 
the eigenvalue problem of each dimer as a non-linear (g)-spin model. 
Analytic investigations of the ensuing three-term non-linear recurrence 
relations are carried out and the respective orthonormal and complete 
eigenvector bases are determined. 

The quantum manifestation of the classical self-trapping in the QDNLS- 
dimer and its absence in the QAL-dimer, is analysed by studying the 
asymptotic attraction and repulsion respectively of the energy lev- 
els versus the strength of non-linearity. Our treatment predicts for 
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the QDNLS-dimer, a phase-transition like behaviour in the rate of 
change of the logarithm of eigenenergy differences, for values of the 
non-linearity parameter near the classical bifurcation point. 
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1 Introduction 



The Discrete Nonlinear Schrddinger and the Ablowitz-Ladik chain models 
have generated intense interest during the last decade to physicists and 
mathematicians. Both systems consist of coupled sets of ordinary differ- 
ential equations, the former has interesting physical motivations while the 
latter has attractive mathematical properties. The property of DNLS that 
makes it physically appealing is that of selftrapping [Jj, ^ |3|, m^.the existence 
of broken symmetry localized states induced by nonlinearity. This property 
makes DNLS an interesting model for discrete solitons or polarons. On the 
other hand the AL equation ||, even though not physically as transparent 
as DNLS has the unique mathematical property of intergrability, as a result 
many of its properties can be investigated analytically 0. 

An issue that already has been addressed by various authors (see e.g. 
15 ||) B 01)' * s that of the quantum properties of these two equations. In 
particular the question of the quantum manifestation of selftrapping and in 
what precise way the quantum versions of DNLS and AL differ. This is 
the question we also address in the present article through the use of novel 
techniques motivated from the area of quantum groups ||. In particular, 
we develop a systematic algebraic scheme for the Quantum DNLS (QDNLS) 
and Quantum AL (QAL) model that uses g-symmetries. This includes the 
finding of the dynamical symmetry algebra for each model and the determi- 
nation of the set of independent constants of the motion for the respective 
Hamiltonians. It is achieved by using a canonical boson realization for the 
su(n) algebra, in the case of the n-site open QDNLS chain. The chain pos- 
sesses then as invariants the set of n — 1 independent Casimir operators of 
the algebra, the eigenvectors of which will determine the linear space of the 
unitary time evolution of the model. These considerations are further ex- 
tended to the case of the QAL chain model. In this case we first show that 
the quantum dynamical variable at each chain site is a particular case of the 
so called g-oscillator || (see [[HJ for related considerations). Then we invoke 
a g-oscillator realization of the g-deformed (quantum) su{n) algebra which 
lead us to determine the invariants of motion of the QAL Hamiltonian that 
are given in terms of the quantum Casimir elements. 

In the second part of the work we turn our attention to the special case 
of two-site (dimer) DNLS and QAL models where the preceding general al- 
gebraic scheme is now applied. In terms of the representation theory of the 
respective dynamical algebras of the two dimers we formulate the eigenvalue- 
eigenvector problem and provide analytic solution in each case. Our aim is 
to use the obtained solutions for a cooperative study of the spectral features 
of the two dimers from the perspective of the phenomenon of seltrapping. 
In particular, we investigate the energy level structure, the precise nature 
of which determines the relative localization-delocalization transition of the 
quantum excitation for values of the non-linearity parameters near the bifur- 
cation of the classical DNLS dimer. The plan of the paper is as follows: In 
the following section we discuss the quantum g-symmetries of QDNLS and 
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QAL, in section 3 we provide the corresponding dimer exact solutions and 
finally in section 4 we give the numerical results with a physical interpreta- 
tion of the resulting plots along with conclusions and future prospects of our 
work. 



2 Quantum symmetries 

2.1 The QDNLS model 

In this chapter we will start the investigation of the dynamical symmetries of 
the quantum DNLS (QDNLS) and the quantum AL (QAL) chain. The two 
models are quantized versions of their classical analogues. The quantization 
procedure is the usual naive method of substituting the classical complex am- 
plitudes by canonical (and non-canonical as in the AL case) boson operators, 
following the normal ordering rule, and the Poisson brackets by commutators 
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We consider first the Hamiltonian of the QDNLS n-site open chain, 



H 



DNLS 



E 



-e(a i+1 + a;_i)c4 - ~(ajaj) 



(1) 



with the canonical commutation relations 



[Ni, cij 
'Ni, a] 



8ij<x\ , 



(2) 



among the boson sites. 

To study the algebraic symmetry of this Hamiltonian we consider the defining 
relations of the su(n) algebra generators, in the Cartan-Chevalley basis \T2 
{ei, fi, , which read 

(3) 
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where r, s are non-negative integers constrained by r+s = 1— a^. Notice that 
the last two equations are the Serre relations which involve cubic products 
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of the n — 1 generators, in terms of these the remaining su{n) generators are 
defined, until the total number of the n 2 — l basis elements is completed. Also 
the Cartan matrix elements which occur as structure constants in eq.(f|), are 



defined as a^- 



25ij 



<5j J+ i — 8ij-i. Finally in this algebra basis the raising 



and lowering generators are conjugate i.e. e\ = /j and the Cartan subalgebra 
generators are self-adjoint viz. h\ = hi, i = 1, ...,n — 1. 

Then we invoke the bosonic realization of the su(n) algebra given in terms 
of the set of n-bosons of eq.(|2f),||l3l, [12| , 



ft 
hi 



a \ a i+l j 



-(Ni - N i+1 ) 



(5) 



Employing the central element h = N1 + N2 + ... + N n (total energy operator), 
we express the number operator at each site in the form 



n-l 



where 



N l = -J2(n-%h l + (n- 1 ) m h, 



1 -1 ... 
1 -1 ... 



n 



1 1 



(6) 



(7) 



is an n x n invertible matrix. 

By virtue of the bosonization of the generators the QDNLS Hamiltonian 
can now be embedded in the su{n) algebra 



H 



DNLS 



E 



-e(ei + fi)-l 



1 n— 1 

z i=i 



ijhj 



x )-h- 



3=1 



2 y 



1\2 



(8) 

This dynamical symmetry in turn implies that all central (Casimir) el- 
ements of the su(n) algebra commute with Hj^^ls-, an d therefore become 
constants of motion. The number of these independent invariants equals the 
rank of the su(n) which is n — 1. From the various forms available of the 



Casimir operators in terms of the su{n) generators ||T3[ we will use one that 
facilitates the comparison with the corresponding g-deformed Casimirs that 
will be discussed next for the case of QAL chain. To this end referring to 
eq.(^|), we introduce the generators 



h a = 



E, 



ab 



z a - e a+1 ei + e 2 H h e n = 0, 

j E ac E cb — E cb E ac a < c < b , no summation 
I e„ a + 1 = b , 



(9) 
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and E a b = 0, for a = b, a > b. Also = E\ h , e\ = e a , and the range 
of values for all indices is from 1 to n, while the index c above is arbitrary 
provided that a < c < b. In terms of these generators the desired Casimir 
invariants are, 

C 2p = (M p ) aa , (10) 

where the matrix M is defined by M = E ac F c b, and recursively the higher 
powers are given be (M p+1 ) a t = (M p ) ac M c b. An alternative set of Casimir 
operators is defined by 

C' 2p = (N^aa (11) 

with N a i, = FacEd). The above expressions provide two sets of Casimir opera- 
tors of any degree for the algebra su(n), i.e. [H DNLS , C 2p ] = [Hdnls, C 2p ] = 
0. The quantum dynamics generated by the Hamiltonian Hdnls is con- 
strained by the existence of these invariances. The state vector of the system 
evolves in the vector space spanned by the eigenvectors of the set of Casimir 
operators that correspond to a given set of eigenvalues. Both the eigenvalues 
and the set of their respective eigenvectors are labeled by a minimal set of 
indices the so called Casimir indices, the number of which equals the rank 
of the dynamical algebra. This is the standard scheme of dynamical symme- 
tries for quantum Hamiltonian systems (see e.g. fl3||). It can be utilized to 
determined the dynamics of the general DNLS chain and shortly it will be 
used in the DNLS dimer case. 

As the spectrum generating algebra of the model, the su(n) algebra, 
provides means for the solution of the eigenvalue problem. This is particularly 
so in the case of two sites to which we now turn. Renaming the generators 
e\ = J + , fi = J_, hi = Jo, in terms of the more usual su(2) spin algebra 
notation, the quantum dimer Hamiltonian after some rescaling and a shift 
by a constant term, reads (for the dynamics of the wavefunction zeros of this 
Hamiltonian see JT3|]), 

H QD = J + + J- + \JI- (12) 

This form of Hq D is particularly suitable for the study of its eigenvalue 
problem which will be taken up in the next chapter. 



2.2 The QAL model 



We turn now to the AL discretization of the continues NLS equation 
This is a Hamiltonian model with Hamiltonian operator 



i 

H AL = -J2{bt(b i +i + b i -i)-2 



i=i 



ln(l + lb\bj 
ln(l + I) 



(13) 



The model employs a set of non-canonical oscillators |J with commutation 
relations 
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b% Oij , 

b\5 i3 . 



(14) 



To demonstrate the quantum group symmetry of the QAL Hamiltonian 
we need first to show the relation of the non-canonical AL-oscillator to the 
so called g-oscillator (quantum Heisenberg algebra) || and subsequently to 
consider the g-bosonization of the quantum group su q {n) . We start with 
the abstract g-oscillator algebra generated by elements a, a', and N with 
commutation relations [PI 



aa' 



qa^a 



q 



-N 



,t 



[N,a] 



-a . 



(15) 
(16) 



The parameter q is the so called deformation parameter and in our case is 
taken to be a non-negative real number. The g-oscillator is so defined that in 
the so called classical limit when q — > 1, it reduces to the standard quantum 
mechanical oscillator. This g-deformed algebra possesses a nontrivial central 
element 



C 



where [N] = (q — q )/{q — q 



-N i 



_1 ), which implies that 



(17) 



a) a = [N] - q N C . 
If we introduce the generators 

&t = <Jq- Nl2 and b = q~ N ' 2 a, 
we express the defining relations in terms of commutators i. e. 



q 



-2N 



N, 6 f l = 6 f , [N, b] = b . 



Then eq. QIq) is equivalent to 

b ] b 

which amounts to 



q- N+1 [N]-gC 



= l- q C 



(1 - q- 2 )tfb . 



(19) 



(20) 



(21) 



(22) 



To correctly identify the AL-oscillator we make the choice C = 0, q 
which imply the relations 



6,6 f ' 



1 + -tfb 
2 



and 



N 



ln(l + Z&t&) 
ln(l + i) 



(23) 
(24) 
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for the g-Heisenberg oscillator. From the representation theory of the ab- 
stract g-oscillator algebra |D|, we know that it possesses a family of irre- 
ducible infinite dimensional representations (irrep), with various operator 
properties assigned to its generators. However the choice C = 0, singles 
out a representation for which the number operator N, is positive definite. 
This Fock-space matrix irrep of the algebra constructed in the Hilbert space 
TCf, spanned by the basis vector {|n >}^L , makes the Hal a proper energy 
operator and provides the following matrix realization for its elements: 



N\n >= n\n > , 
b\n >= \J{n}\n - 1 > 



(25) 
(26) 

&t| n >= ^{n + l}|n+ 1 > , (27) 

with 6|0 >= and {n} = £f . 

We proceed now showing the quantum group symmetry of the QAL- 
chain. First let us recall the defining relations of the quantum su q (n) alge- 
bra according to the Jimbo-Drin'feld scheme, in the Cartan-Chevalley basis 

[e h f h h = q hi Y i=1 , they are § 



ki 
k^ kj 

[e;, fi) 



k% k{ — 1 
hk- 

I-.. 

[2/cj]5y , 



(28) 
(29) 



together with the g-Serre relations 

1 — a 



E(-i: 

r,s 

E(-i)' 



1 O^ij 



C 2 C J C J 



J q 



£T £ -fS 



J g 








i 7^ j- 



(30) 



The last ones involve the so called g-binomial coefficient [[16] defined in terms 
of g-factoral, [m]\ = [l][2]...[m] by 



m 
n 



[m]\ 



[n] \[m — n] 



(31) 



Let us mention that here we consider the su q (n), only at the algebra level, 
the rest of its Hopf algebra structure is not needed for our present purposes, 
but they remain worth studying from the physical point of view. By analogy 
with the bosonization of the su(n) by a set of canonical oscillators, the q- 
bosonization of the quantum su{n) is always possible by using the g-deformed 
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oscillator of eq.([L6|). Indeed for a set of n such oscillators the bosonized 
generators are |T7[ 

ei = a\di + i , 
f — t 

jfe. = gW--^!). (32) 

We also introduce the generators C{ = q 1 / 2 ( Ni+Ni + 1 \ % — 1 ( 2, n — 1, which 
commute with any {e^, /j, fcj} with the same index. Using this commutativity 
property and eq.fll9f) we express the su q (n) generators in terms of the QAL- 
bosons i.e. 

fi = b\ +l biCiq- 1 , 

By means of these equations and the relations (|24j ) we now write the model's 
Hamiltonian as 

n 

Hal = ~Y, {qCr\ei + fi) - 2N,} , (33) 

i=l 

Due to eq.(||), we can finally embed this Hamiltonian in the su q (n) algebra 
as 

n [ n— 1 ] 

H al = -E <l C i\ & i + fi) ~ E(^%^ - (V-^inh (34) 

i=l [ 3=1 J 

where 

— y — q ^ ' ■ 

(35) 

It is important at this point to emphasize that as our analysis shows 
that the quantum group symmetry of the QAL model stems from the non- 
canonical character of the quantum variables defined at each site, (c.f. eq.(|i~4])). 
This non-canonical character of the degrees of freedom is in fact necessary in 
order to prove the Hamiltonian structure of the AL equations with respect 
to the Hamiltonian of eq.([L3|), ||. We conclude therefore that the found 



quantum group symmetry is a genuine feature of the model 

Similarly to the preceding case of the DNLS su(n) invariances, the su q {n) 
dynamical symmetry of the QAL chain model implies that all quantum 
Casimir elements of su q {n) commute with Hqal- In this way they become 
constants of the motion generated by the QAL Hamiltonian. For the con- 



struction of g-invariants of su q (n) there are several approaches 



For our purpose a construction similar to that outline before for the DNLS 
case will be employed pi |. Let us introduce in terms of the g-deformed 



algebra generators the following elements: 

h a = e a - e a+ i ei + e 2 H h e n = 0, 

E ac E c b — q' 1 E C0 E a c a < c < b , no summation 
E ab = { e a a+l = b, (36) 

(q-q' 1 ) a = b, 



7 



and E ab = for a > b. Also F ba = El b , e\ — e a , and the range of values for 
all indices is as before. 

Notice that the g-deformed generators E ab (q), are not invariant under q — ► 
g -1 , so we can define another set of generators E ab = E ab (q~ l ), F ab = 
F a b{q 1 )- Then we introduce the elements 

M ab = E ab q^- 2a F cb q^ +2c , (37) 

which yields the g-Casimir elements 

C 2p = q 2a (MV) aa . (38) 

By the exchange q — > g -1 , we obtain a new set of invariants 

C 2p = C 2p {q- 1 ) = q- 2a (MP) aa , (39) 

where M p+1 as before is obtained recursively from the operator- valued matrix 

M ab = M ab (q- 1 ) = E ac q-^ +2a F ab q-^- 2c . (40) 

Finally two Casimir generators invariants under the exchange q — ► q~ 1 
may be constructed by combining the two previous ones in the form 

C 2p = {C 2p + C 2p )/{q + q- 1 ) 
C 2p+1 = (C^-C^/iq-q- 1 ). (41) 

These last expressions of the invariant Casimir operators modified with ap- 
propriate coefficients reduce presicely in the q — > 1 limit to the Casimir 
operators of equal degree for the non-deformed su(n) algebra (see [21], for 
details and some additional expressions for g-Casimir elements). 

If we now confine ourselves to a two-site restriction of the QAL-chain 
then the Hamiltonian 

Hal = -qC^(a\a 2 + ajoi) + + N 2 ), (42) 

has a su q (2) dynamical symmetry. This symmetry will be utilized in the next 
section in order to solve the eigenvalue problem of the AL-dimer. 



3 Exact solutions 

3.1 The DNLS dimer 

We now turn to the analytic determination of the eigenvalues and eigenvec- 
tors of the quantum DNLS chain in the simplest case of two sites. As was 
already mentioned there is a number of previous investigations in the case of 
classical dimer, which predict the well known phenomenon of selftrapping Jl| . 
But also for the quantum case there are numerical and perturbative efforts 
0, and quasiclassical treatments along with formal techniques based on 
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the inverse scattering method ||. aiming to study the quantum manifesta- 
tions of selftrapping, 

Our approach will be entirely analytic and will be based in the su(2) 
pseudo-spin expression given to Hqd Hamiltonian (0). According to pre- 
viously given general prescription of the group symmetry of the chain, the 
two-boson realization of the su{2) generators reads 



J+ — 01^2, 

Jo = -(iVi-iVs). (43) 

The generators obey the standard relations of the quantum angular momen- 
tum 

[JoiJ±\ = =t^±j 

[J + ,J_] = 2J . (44) 

With the Casimir operator adapted from the preceding general formulae 

C = J (J - 1) + J+J- = Jo{Jo + 1) + J-J+ , (45) 

the matrix form of the Hamiltonian is induced by the (2j + l)-dimensional 
matrix representation of the generators in the spin basis {\jm >Y m= - where 

3 = 21 ^-i I'" - ' rea( is: 

C\jm> = j(j + l)\jm>, 

Jo\jm > = m\jm >, 

J±\jm > = J (j =F m) (j ± m + 1) \jm± > . (46) 
In the two boson labeling of the spin state vector 

\j = -(ni + n 2 ) ; m = -(m -n 2 ) >= |ni,n 2 >, n 1(2 > 0, (47) 

we see that the total boson excitation number n\ + n 2 , determines the di- 
mensionality of the Hamiltonian matrix, and also the total energy allocated 
to the dimer. Then the eigenvalue equation 

H QD \ ( f > >=\\4 > > (48) 

for a Hamiltonian eigenvector which is expressed by 

j 

|0 >= °n\ n >> ( 49 ) 

n=-j 

leads to the difference equation 

7n 2 c n + e (j + n)(j - n + l)c n _i + \J (j - n)(j + n + l)cw+ij = Ac n . 

(50) 
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with boundary conditions c_j_i = Cj + % = 0. We introduce new variables 
"071 = c„7 n and impose upon 7„'s the relation 



7n\/ '{] + n)(j - n + 1) = 7„_i (51) 

with solution 



J-j+n 



\ 



n-l 

n!n(2j-*)T-i, n=l,2,..,2j. (52) 

i=0 

In this way we obtain a normalized form of the iteration 

ipi-j+i = (A + + /ifVi-j-i , i = 0, 1, j, (53) 

where 

A4 1} =7(i-i) a (54) 

and 

fj£ 2) = -i(2j + l-i). (55) 

(also we set /4 = 1 for later use). Then from the boundary conditions 
ip-j-i = ipj+i = and the seed of the recurrence relation ip_j = 1 , we 
deduce that the polynomial p(X) = ipj + i(X), is the characteristic polynomial 
of the Hamiltonian matrix Hqo- Let us write 

P(A)= Ef^ 2 - (56) 

i=0 



then the coefficients of the characteristic polynomial are determined by in- 
duction. Ii 
defined as 



duction. Indeed if we introduce the concept of the weight of a coefficient 5- 



wt(5? } ) =Wi = degp(X) - i = 2j + 1 - i, (57) 
then we shall obtain the following expression for the coefficients 

W - EE E 

ii, 12,... ,1^=0 Z Wi =ii + ... 4-itoj— 1 iu) i _ 1 =ii+..-+iui i _ 1 — 1 
'3— «3 '2— «2 . . . 

e e /^-MirV!; , (58) 

(2=11+12—1 il=*l — 1 

where the summation indices satisfy the constrains 

i w . ^ 2, ii + i 2 + ...i^ = (59) 

and if ifc = 2 then z^+i = 0. 

Next we determine the associated eigenvectors by specifying their expan- 
sion coefficients in the spin basis (c.f. eq. (f49|) ) . Since the n-th of these co- 
efficients is identified with the n-th order iteration of the recurrence relation 
issued in eq.(|53"D we write 

4>-j+n(Xa) = q{K)- (60) 
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The polynomial q(X a ) for each of the eigenvalues {A o }^_ , has degg(A a ) = n, 
then we can write 

n 

9(A a ) = E4 n " 1) A^. (61) 

The finding of the eigenvector amounts then to the determination of the 
coefficients 5^ . This is done by use of eq.([53|) and the solution is again 
given by eq. (|58|) , with = 1. The variable u>j this time stands for the 
weight of the coefficients of the q(X) polynomial and it is defined by 

wt(5f n ^) = Wi = n — i , i — 0, 1, n . (62) 

Finally, we shall show that the found operator spectrum possesses proper 
mathematical properties, namely that the energy eigenvalues are real num- 
bers with zero multiplicity (non-degenerate) and that the corresponding set 
of eigenvectors form an orthonormal and complete basis. To this end we 
first recall some generalities ||22|| . Let us assume that we have a set of single 
variable polynomials {q n (x)}^~Q, obeying the three term recurrence relation 

a n q n+1 (x) + (3 n q n (x) + 7„(x) = xq n (x) , (63) 

with g_i(x) = 0, which have roots x iy i.e. p n (xi) = 0, if we introduce the 
the numbers d n = J-^^d^i, with do arbitrary positive number and the 

variables a n = a n a n+1 , then the socalled Darboux-Christoffel (DF) formula 
given by 

y- 1 qn(x)q n (y) = On-i 1 q N (x)q N -i(y) - q N ^i(x)q N (y) 
d l a N 4r-i x-y 

is valid for the polynomials. For two roots x = Xi,y = Xj the DF formula 
becomes 

" W Qn(x)q n (y) 



E = M? (65) 

M i = E — n~ = * — QNWqN-iixi) ■ (66) 

n=0 a n a N «JV-1 



with 



The non-degeneracy of the roots is now obtained from eq . (|65|) , which implies 
that for Xi = Xj, we get q N (xi) ^ 0. Turning now to our case we can write 
for the H DNLS eigenvectors the expression 

m\a) >= K E t±^L\ - J + k>, (67) 

k=0 6k 

where the normalization constant (which determines the value of j-j) 

— 72 — > ( 68 ) 

- /V a k=0 € k 
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involves the factor 

' 2j 



(69) 



with a = 0, 1, ... , 2j. 

To verify the orthonormality of the eigenvectors we resort to the quoted DF 
formula and write 

m\ a m\ b )) =Af a Kj: ^+*(W-*-+*( A ») =Sab (70) 

k=0 € k 

We note that the orthonormality of the eigenvectors is equivalent to the 
orthonormality of the coefficients of the eigenvectors in the spin basis taken 
as polynomials with discrete variable, which is the corresponding eigenvalue. 
Also the completeness of the eigenvalue basis can be based on the DF formula; 
indeed it is straightforward to write: 

E Wa.))Ma.)| = E (E^ ^^f^^ l H+k)(-j + i\- 

a=0 k,l=0 {a=0 6 k 6 l J 

(71) 

Having constructed the eigenvector basis of the Hamiltonian, we can now 
investigate dynamical questions, by expanding a suitable initial state vector 
on this basis, here however we refrain from doing so and shortly we will turn 
to the study of the quantum manifestations of the classical selftrapping. 



3.2 The QAL dimer. 

In this section we put forward an analysis of the eigenvalue problem for 
the QAL-dimer based, as in the DNLS-dimer, on the symmetry algebra of 
the Hqal Hamiltonian. The new feature in this case is that we deal with 
the quantum, or g-deformed algebra su q {2). As the value of g-deformed 
parameter is a non-negative real number, the representation theory (which 
as we have shown in the preceding case is an indispensable element for the 
concise formulation of the eigenvalue problem) of su q {2) is very similar to its 
q — 1 limiting case. 

Starting with the Hamiltonian of eq.p2]) and after appropriate shifting 
and scaling with constant factors we get 



Hal = E a l( Q i+i + a i-i)- ( 72 ) 
=i 



This can be expressed in terms of the effective quantum su{2) generator ( 
viz. 

H AL = J q + + Jl. (73) 

(We note here that the Hal Hamiltonian is exactly that of Azbel-Hofstadter 
model which describes Bloch electrons in magnetic field p3 , however in that 



case the g-deformation parameter is determined by the value of the magnetic 
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flux and is root of unity; given that the latter is solvable by Bethe ansatz 
techniques the relation between the two models is worth studying.) These 
quantum generators bosonized in an angular momentum fashion become 

J% = dial, 

4 = l(N!-N 2 ). (74) 

They satisfy the quantum su q {2) algebra commutation relations i.e. (q = e s 

), 

[JlJl] = ±4 

{JlJl] = [2 Jq9] = su^l) (75) 

With the g-Casimir operator adapted from the preceding general formulae 
for g-deformed central elements, 

C q = [4] [J Q o ~l] + JU q - = [4] [4 + 1) + J1JI, (76) 
the matrix representation of the q-deformed generators 
C q \jm> = + l]\jm >, 



Jl\j,m> = ^m)\j ^m + l]\j,m+ 1 >, 

Jo\h m > = m\j,m>, (77) 

is very similar to their q = 1 counterparts. We notice that the non-linearity 
of the model resites in the fact that the off-diagonal matrix elements of 
the Hamiltonian matrix are expressed in terms of q- numbers (c.f. eq. (|77|) ) . 
Simple inspection shows that in the limit 7 — ► 0, q — > 1, where for any q- 
number the limit [x] — > x, is valid, the generator matrices J± become the 
su(2) generators i.e. J± J±. In this limit the QAL-dimer becomes a linear 
angular momentum problem that is easily solved. 

To address the eigenvalue problem we first note that a similar relation 
as in eq.(f4~?D for su(2), holds between the state vectors of su q {2) and their 
parametrization in terms of the excitation numbers of two g-bosons. If we 
now assume an expansion \<f> >= J2 J n =-jCn\n >, for the eigenvector \<j) >, 
associated to the eigenvalues A i.e. Hqal\(P >= >> then this last relation 
by means of eqs.flIB[),(|77D yields 



\J[] - n + l][j + njc n _i + \J[j + n + l][j - n]c n _i = Ac n . (78) 

Again we normalize this difference equation with a change of variables c n = 
— where the 7 n 's satisfy the relation 



7n 



7n = 7 n -iy/\j + n]\j-n + l], (79) 
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and explicitly are given by 



1-j+r. 



\ 



n-1 



H!ll7-j > n=l,2, ... ,2j. (80) 



i=0 



Then the normalized iteration with vanishing boundaries z. e. Vj+i = V'-i-i = 
and seed V'-j = 1 reads, 

VWi = - [j + n][j -n + l]V> n -i- (81) 

As was explained in the similar situation of the DNLS-dimer, the eigenvalues 
are the roots of the real polynomial ipj+iW — 0. By induction we find 
that one needs to distinguish between integer and half-integer values for the 
index j. In other words we need to distinguish between even total number 
of excitation quanta and odd total number of excitation quanta that are 
available in the two g-boson sites of the model, (c.f. eq. (f£7|) . Explicitly we 
find that for j-integer 

^• +1 (A) = Ag(A 2 ) = (82) 

where degg(A) = j and degV^+i (A) = 2j + 1. This implies that the energy 
spectrum contains the zero eigenvalue and that the remaining energy levels 
are arranged in doublets for each pair ±A. If we now define x = A 2 , and 
assume that 

Q(X) = t ^ J) X m , (83) 

m=0 

then we seek to determine the coefficients (3^ J \ 
Similarly for j-half-integer we obtain that 

^■ +1 (A)=p(A 2 ) = 0, (84) 

where degp(A) = 2 4^-, and 

|(2i+i) 

p(*) = E PSPxT- (85) 

m=0 

Then by induction of the recurrence relation of eq. fl81~|) , we specify the (j + 1)- 
order of iteration ^j +1 (A), which provides the coefficients 

v ij— m — 2 23 — 2 i2 — 2 

Pm = E E E E fki-m-VtoHil- (86) 

ij-m=Pj-m ij-m-l=Pj-m-l «2=P2 »X=Pl 

In this expression we have used the abbreviation 

fin = -{] + n]{j - n + I], (87) 
and the p r 's are indices determined by the iteration 

p r = Pr ^ + 2 r = 2,3,... , Pl = l. (88) 
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Next we construct the set of (2j + l)-eigenvectors of Hqal by specifying their 
expansion coefficients ip_j +n , n = 0,1, 2j, in the g-spin basis. We find that 
for n-even 

4>- j+n (\) =u(X 2 ), (89) 



and 

Similarly for n-odd 
with 



u 



w 



(x) = E ^- 1] x m - 

m=0 

^_ i+n (A) = w(X 2 ), 

l/2(n-l) 

(x)= E P%- 1] x m 

m=0 



(90) 
(91) 
(92) 



where the coefficients S^~^ are determined again by the relation (]86f) . 

Concerning normalization of the eigenvectors we note that since we have 
set ip-j = 1, the remaining unspecified variable 7_j, is determined from the 
normalization condition of the eigenvectors as in the preceding case. The 
latter are given explicitly by 



tp- j+k (\ a ) 



k=0 Ek 



j + k>, 



(93) 



where the normalization constant 

1 



27 f-^W,) 

k=0 



involves the factor 



for a = 0,1, 2j. 



efc = S 0)k +[k}\, 



k 



(94) 



(95) 



J g 



4 Numerical Results and Discussion 

We use now the previously derived exact results to study the spectral be- 
havior of QDNLS and QAL models. These results are shown in Figure 1 
for QAL and in Figure 2 for QDNLS respectively. In Figure la we present 
the spectrum of QAL as a function of nonlinearity for j = 4. We note the 
distinct "repulsion" of the levels: the larger the value of nonlinearity the 
greater the repulsion between the levels. This tendency for repulsion of the 
energy levels is not related to similar repulsive and avoiding crossing behavior 
one encounters in quantum non-integrable systems (see e.g. f24|), since our 
QAL-dimer is an integrable model. It should rather be attributed to the non- 
linear dependence of the eigenvalues on the g-deformation parameter, that 
is to say that it should be an effect of the quantum group symmetry itself. 
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In Figure lb we plot the lowest QAL energy levels for a given nonlinearity 
value (7 = 2) but as a function of the number of quanta in the system. The 
classical AL is obtained in the limit of very large j-values. We note here the 
repulsive aspect of the levels as well. 

In Figure 2a the spectrum of QDNLS is shown for j = 3 as a function 
of the DNLS nonlinearity parameter. We note a great difference from the 
corresponding spectrum of QAL; here pairs of adjacent levels are grouped 
together and they merge as nonlinearity increases. This property is evident 
for all pairs of levels and especially in the lowest ones. We note that the 
rate of convergence is 7-dependent. In order to investigate this property we 
plot in Figure 2b the energy difference of the lowest two pairs of levels as 
a function of the nonlinearity parameter in a log-log plot. The labels (1) 
and (2) signify the lowest and next to lowest pairs respectively. There is a 
drastic change in the convergence plot as the amount of nonlinearity exceeds 
a certain value. While at small 7 the energy difference shrinks very slowly, 
the rate changes very rapidly for large 7. The convergence is done seemingly 
in an algebraic fashion with energy dependent exponents that can be easily 
obtained. 

The numerical results obtained for the energy splittings demonstrate the 
role of classical selftrapping in the quantum mechanical regime. In its ab- 
sence, as in the AL equation, the system energy spectrum shows repulsion, 
while in its presence the opposite effect is manifested, viz.level clustering. 
This energy clustering in the quantum regime is a signature of long-lived 
excitations. Indeed, if an excitation is created on one site of the QDNLS 
dimer, while the latter is in a large nonlinearity regime, the tunneling time 
will be very large, leading thus to a very long lived and localized excitation. 
This excitation can be thought of as a quantum breather induced by the 
nonlinearity of the corresponding classical problem. Such a breather does 
not exist in the QAL becomes evident from its spectrum. Both clas- 

sical counterparts of the two problems studied in this work are integrable; 
it would be interesting to investigate in detail and compare cases with and 
without selftrapping but while the system is also non integrable. In such 
problems the simultaneous presence of classical localization and chaos will 
have interesting reprecautions in the quantum mechanical spectrum. 
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Figure captions 

Fig. la. The spectrum of QAL as a function of nonlinearity for j = 4, i.e. 
nine energy quanta. 

Fig. lb. The lowest QAL energy levels for a given nonlinearity value (7 = 2) 
as function of the number of quanta in the system. 

Fig. 2a. The spectrum of QDNLS is shown for j = 3 as a function of the 
DNLS nonlinearity parameter. 

Fig. 2b. The difference of the lowest two pairs of energy levels as a function 
of the nonlinearity parameter in a log-log plot. The labels (1) and (2) signify 
the lowest and next to lowest pairs respectively. 
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